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1.  Introduction 


12  3 

Since  Marsh  '  '  proposed  his  acoustic  scattering  concept, 

4 

the  concept  had  been  generally  misunderstood.  This  author  had 
extended  his  theory  to  imperfect  reflectors  and  found  the  theory 
very  useful  and  physically  revealing. 

5  G 

Recent  translations  '  from  Russian  literature  show  that  the 
Russian  community  has  analyzed  acoustic  wave  scattering  by 
utilizing  Cass's  perturbation  method.  The  method  of  Bass  was 
initiated  in  19C0  for  the  community  of  electromagnetic  field 
about  the  time  Marsh  formulated  his  for  the  community  of 
underwater  acoustic  field.  The  Russian  community  has  since  been 
aware  of  the  concept  of  Marsh-Kuo. 

Somehow,  Bass's  method  got  quite  popular  in  the 
England.  It  became  the  basis  of  some  publications^, 
reporting  solutions  to  problems  basically  addressed 
of  Marsh-Kuo  ten  years  earlier. 

In  this  letter,  an  attempt  is  made  to  identify  the  common 
origin  of  these  scattering  concepts. 

2.  A  Genetic  Tree  of  Green 


U.S.  and 
9  10 

by  the  method 


Let  us  start  from  the  Green's  identity  (see  e.g.  Baker  and 
Copson^^)  which  takes  the  following  form: 

JJJ  Jj <is' 

V  s' 

where  S'  is  a  closed  surface  bounding  the  volume  V  and  &/9N 
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denotes  differentiation  along  the  outward  normal  to  S ' .  This 
identity  is  valid  when^and^and  their  first-  and  second-order 
partial  derivatives  are  continuous  within  and  on  S'.  In 
particular  i f  ^and ^ sat  is f y  following  Helmholtz  equations: 

-C)  -  -4*rr©.c-?i'> 


Where  =  W/c  ,  x  =  (x,  y,  z)  ,  £'  =  (x ' ,  y',  z') 


then,  equation  (1)  is  reduced  to  the  Green's  formula  (e.g.  Bass 

5 

and  Fuks  )  : 

<<>«>•  J[(  ft 


In  the  first  integral  the  region  of  integration  is  the  Volume  V, 
filled  by  sound  sources.  In  the  second  integral  the  region  of 
integration  is  the  surface  S',  filled  by  a  certain  distribution 
of  simple  and  double  sources  represented  by  4^  and 
respectively. 


When  C  =  0,  equation  (4)  becomes: 


S' 


(5) 


in  wliich  is  the  location  vector  for  the  bounding  surface 
points.  The  bounding  surface  S'  can  be  divided  into  a  local 
boundary  surface,  s,  and  a  hemi-spher ical  surface  at  infinity 
(say  ).  When  Sommerfield  conditions  are  satisfied  (see  e.g. 
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Baker  and  Copson^^),  equation  (5)  is  reduced  to  the  following 
Green's  formula: 

To  summarize,  the  genesis  of  equation  (6)  is  the  Green's 
identity.  Equation  (6)  is  a  specialized  form  of  the  Green's 
identity  for  a  wave  field  when  sources  are  distributed  only  on  a 
surface  s.  It  simply  says  that  the  velocity  potential  field  at 
any  point  in  space  is  defined  by  the  potential  value  and  its 
normal  derivative  on  the  surface  s.  Equation  (6),  in 

turn,  can  be  identified  to  be  the  genesis  of  existing  surface 
scattering  theories,  while  equation  (4),  without  surface  sources, 
is  the  genesis  of  existing  volumetric  scattering  theories. 
Starting  from  Green's  identity  through  Green's  formula  to  wave 
scattering  theories,  it  forms  a  kind  of  a  genetic  tree,  say 
Green's  Tree,  of  the  origin  of  various  scattering  concepts. 

3 .  Identification  of  Past  Acoustic  Scattering  Concepts 

Based  on  Green's  formula  (6),  scattering  concepts  for  rough 
boundaries  can  be  identified  with  modern  existing  theories. 

Let's  identify  the  geometry  of  interest  as  in  Figure  1.  A 
general  class  of  Green's  function  can  be  given  by  (e.g.,  Morse 
and  Feshback^^): 

^  9.cz,s‘> 


Where  =  (x  -  x' )^  +  (y  -  y*)^  +  (z  -  z*)^. 

According  to  properties  of  Green's  function,  the  first  term 
is  the  proper  form  of  Green's  function  for  a  three  dimensional 
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point  source  in  free  space.  The  second  term  accounts  for  the 
effect  of  boundaries. 

Given  an  incident  wave  field  in  z  >  z*  ,  a  scattering  tlieory 
seeks  to  estimate  the  effects  of  a  rough  boundary:  =’J(x*  ,  y*  ). 

The  incident  wave  field  is  assunicd  effectively  plane,  thus  its 
effective  source  location  is  at  oo  .  The  Green's  formula, 
equation  (6),  asserts  that  the  reflected  wave  field  can  be 
obtained  by  summing  the  effective  surface  sources  induced  by  the 
incident  wave  field  on  the  rough  surface  s,  according  to  the 
integral  formula.  In  that  integrand,  ^  represent  the  total 
field  in  the  upper  medium  (z  >  z')  evaluated  at  the  rough 
boundary,  while  Green's  function  ^  “  V*  ^ 

solution  to  the  differential  equation  (3)  and  is  utilized  to 
express  the  directivity  pattern  of  these  effective  surface 
sources . 

In  a  similiir  way,  another  Green's  integral  formula  for  the 
transmitted  wave  field  in  the  lower  medium  (z  <  z')  can  be 
derived.  It  will  be  of  the  same  form  as  equation  (C).  However,  ^ 
in  the  integran<l  therein  will  be  the  transmitted  wave  field 
evaluated  at  the  rough  boundary  and  Green's  function  will  satisfy 
a  differential  equation  like  equation  (3)  with  appropriate  k  for 
the  lower  medium. 

Basically  then,  the  above  discussion  implies  that  the  rough 
boundary  can  bo  replaced  by  a  distribution  of  surface  sources  of 
the  appropriate  strength  and  directivity.  In  this  respect,  these 
effective  surface  sources  are  in  free  space  and  the  form  of 
Green's  function  without  in  equation  (7)  is  quite  general  for 
the  study  of  acoustic  wave  scattering  at  rough  boundaries. 

By  assuming  Green's  function  in  the  spectral  forms 


JIf  e 


(B) 


Where 


Hi,, 

A  (K,  x')  is  found  from  equation  (3)  to  be; 


e 


(9) 


i  .  e . 


X-  Uf  *£ 

J)J 


(10) 


It  can  be  shown  that  integration  of  the  above  equation 
produces  the  result: 


e 


la-*'  I 


(11) 


This  Green’s  function  in  three  dimensional  free  space  can  also,  be 

obtained  from  tlie  physical  nature  of  Green’s  function  (see  e.g. 

1 2 

Morse  and  Feshback  ).  The  method  selected  liere  is  to  illustrate 
the  tie  between  modern  scattering  theory  and  Green's  formula. 

It  can  also  be  shown  that  integration  with  respect  to  k  only 
yields  the  result: 
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(12) 


for  z  -  z*  >  o  cjnd  z  -  z^  <  o,  respectively.  That  is 


(13) 


for  z  -  z^  ^  O  and  z  -  z*^  <  O,  respectively.  A  different  method 

13 

was  utilized  by  Drekhovskikh  to  derive  a  similar  equation  for 
which  z^  =  o. 


Noting  that  the  inward  normal  n  (=  -  &)  in  Figure  1,  is  given 


by 


(14) 


the  following  expression  can  be  obtained  from  equation  (13); 

(»-a'i  “  *it  JJ  ® 


D 


-  rK*'-*  W.'l'.*'  3  — 


(15) 


for  z  -  z^>0,  i.e.  the  upper  medium. 


The  following  relationship  transforms  ds  to 


dA  =  dx  ^  dy  *  : 
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(16) 


Utilizing  equation  (11),  the  Green's  formula  (6)  can  be 
written  as: 


_ _ 

I 


(17) 


This  is  the  well  known  theorem  of  Helmholtz  or  the  Helmholtz 
integral  (e.g.,  Baker  and  Copson^^).  Tliis  integral  form  of 
solution  has  been  a  popular  starting  point  of  published  studies 
in  high  freque?ncy  scattering  by  using  Kirchhoff's  local  plane 
approximation  (see  e.g.,  Eckart^^,  and  Mintzer^^). 


Utilizing  equation  (13),  (15),  and  (10),  equation  (17)  takes 

the  following  form: 


4^fs>=  J[  e 


F(^j  )  d 


(IB) 


V^here 
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(19) 


e 


J»5 


By  defining  the  following  Fourier  Integral  Transform 


/w  *  I  rr  *•‘^’‘^^'4'^') 

^iry-  JJ  ^ 

»>V 


(20) 
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equation  (16)  cein  be  written  as 


(21) 


It  is  to  be  noted  that  F  (k  ,  k  )  is  a  Fourier  Integral  Transform 

X  y 

of  a  complicated  statistical  function  of  the  roughness 
character  i  Htica  roproaanted  by  quantitloa  like  ,  and 

of  the  boundary  properties  indirectly  represented  by  quantities 
like  ^  {y*  ),  ( jc  *  )/3  x*  ,  that  are  affected  by  the  densities 

and  acoustic  velocities  of  the  contacting  media. 


According  to  equation  (13)  representing  a  quasi-plane  wave 
expansion,  equation  (21)  implies  the  following.  The  total 
contribution  from  an  equivalent  plane  source  distribution  given 
by  F,  equation  (19),  to  the  reflected  wavefield  in  the  upper 
medium,  (z  >  z*  )  at  an  observation  point  can  be  expressed  as 
the  combined  effects  from  quasi-plane  waves.  Amplitudes  of  these 
quasi-plane  waves,  F,  are  a  function  of  radiation  directions  (k^, 
ky),  boundary  characteristics,  and  the  incident  wave  field.  They 
are  quasi-plane  wave  because  many  radiation  directions  are  of  the 
imaginary  decaying  type. 


A  similar  process  will  lead  to  an  integral  form  similar  to 
equation  (21)  for  the  transmitted  wave  field  in  the  lower  medium 
(z  <  z'  ) . 


These  are  exactly  the  plane  weive  expansion  forms  utilized  as 
a  starting  point  for  the  development  of  surface  scattering  theory 
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by  Marsh  et  al  '  '  and  Kuo  .  Thus,  equation  (21)  is  the  Quasi 
Fourier  equivalent  of  the  Helmholtz  integral. 

In  the  form  of  equation  (21),  the  effects  of  incident  wave 
and  boundary  characteristics  (both  roughness  and  acoustic 
properties)  on  the  reflected  field  are  seen  to  be  expressed  in 
terms  of  their  total  effects  on  each  Generalized  Fourier 
Component  of  the  reflected  field.  Each  Generalized  Fourier 
component  represents  a  reflected  field  in  a  different  reflected 
direction.  At  this  point,  the  perturbation  concept  of  Marsh-Kuo 
proposed  the  expansion  of  this  random  Fourier  component  amplitude 
in  terms  of  the  pertinent  statistical  properties  consistent  with 
boundary  conditions  of  interest. 

Another  basis  to  modern  surface  scattering  theories  came  from 

7  5 

Bass  '  and  can  be  divided  into  two  types. 

7 

The  original  scattering  concept  advanced  by  Bass  was  to 
express  the  reflected  field  as  the  sum  of  mean  and  stochastic 
fields.  Then  he  applied  small  jier  t  urba  t  ion  to  and  performed 
statistical  average  of  the  appropriate  boundary  conditions.  The 
latter  process  produced  two  boundary  equations  or  conditions. 

The  first  equation  was  the  mean  field  boundary  equation  in  which 
the  mean  field  at  the  boundary  was  related  to  mean  statistical 
properties  of  the  boundary  and  to  correlation  between  random 
roughness  and  the  stochastic  field  at  the  boundary.  The  second 
equation  was  the  stochastic  field  boundary  equation  in  which  the 
stochastic  field,  at  the  boundary  was  related  to  random  boundary 
roughness  variables  and  the  mean  field  at  the  boundary.  This  - 
stochastic  field  at  the  boundary  was  identified  as  the  surface 
boundary  source  to  the  stochastic  field  in  space  through  the  use 
of  the  Green's  formula  (C)  with  an  appropriate  Green's  function 
evaluated  at  the  mean  surface.  The  resulting  Green's  integral 
equation,  then,  expressed  the  stocahstic  field  in  space  in  terms 
of  a  mean  field  at  the  boundary  and  random  boundary 
characteristics.  By  utilizing  this  stochastic  integral  equation. 
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the  stochastic  field  in  the  mean  field  boundary  equation  (the 
first  boundary  equation)  can  be  eliminated.  Thus,  the  mean  field 
at  mean  boundary  is  an  integral  representation  like  the  Helmholtz 
integral  whose  surface  source  is  a  function  of  the  mean  field  and 
mean  roughness  properties  at  the  mean  boundary.  Based  on  this 
Helmholtz  integral  representation  of  the  mean  field  and 
identifying  the  mean  field  as  that  of  the  well  known  field 
produced  by  the  plane  boundary,  quantities  such  as  reflectivity 
can  be  computed. 

The  above  type  of  concept  formed  the  basis  of  problems  of 
perfect  reflectors  reported  by  Wenzel^  and  by  Dowling  et  al^^. 

Later  in  his  book^,  Cass  reported  a  change  in  his  method 

(possibly  reported  originally  by  Alekhin^^).  The  only  change  was 

to  use  the  plane  wave  exjiansion  form  as  (21)  to  continue  his 

stochastic  field  at  the  boundary  into  space.  Thus  the  new  method 

incorporated  the  Quasi-Fourier  type  approach,  or  the  plane  wave 

expansion  approach.  Analyses  based  on  this  type  of  approach  were 

17  9 

publishe?d  by  Lysanov  and  Kuperman  . 

The  modified  Bass  method  or  even  ttie  old  method,  is  not 

really  different  from  that  of  Marsh-Kuo  in  perturbation 

principle.  The  latter  method  starts  from  expressions  like  (21) 

for  reflected  and  transmitted  fields.  Then  F  (k  ,  k  )  for  each 

X  y 

fields,  was  found  by  formally  satisfying  the  required  boundary 

conditions  to  the  second  order  in  the  boundary  roughness 

characteristics.  F  (k  ,  k  )  for  each  field  thus  found  is 

X  y 

substituted  back  into  the  form  of  solutions  like  (21).  The 
results  are  the  total  reflected  and  trnsmitted  fields.  A 
statistical  average  of  those  results  may  be  identified  with  the 
averaged  field  solution  of  Bass  et  al. 

So  far,  tlie  assumed  form  of  solution  (21)  utilized  by 
tlarsh-Kuo,  was  shown  to  be  the  Quas i -Four i er  equivalent  of 
Helmholtz  integral  and  as  such  an  exact  solution.  It  has  been 
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also  briefly  indicated  how  the  Bass  scattering  theories  are 
genetically  located  in  Green's  tree  relative  to  Marsh-Kuo's 
theory.  In  principle,  one  should  be  able  to  show  how  they  are 
related.  However,  because  of  different  approaches,  fully 
spectral  approach  of  Marsh-Kuo  on  one  hand  and  fully  spatial  or 
mixed  spatial-spectral  approaches  on  the  other  hand,  assessment 
of  the  equivalence  is  difficult  without  detail  analysis.  This  is 
not  in  the  scope  of  this  letter.  However,  we  shall  show  in  the 
following  the  equivalence  to  the  first  order  in  the  boundary 
characteristics. 

Dowling  et  al^^  considered  an  incident  wave  of  time 
dependence  exp(-iA^t).  It  was  assumed  to  propagate  in  the  general 
direction  against  vertical  coordinate  z  and  in  the  same  general 
direction  of  horizontal  coordinates  (x,y).  The  scattering 
phenomenuni  of  this  wave  from  an  underlining  rigid  and  rough 
boundary  located  at  z*0,  was  investigated. 

By  following  Bass,  the  total  potential  ^  was  divided  into 
a  mean  field  ^  and  a  stochastic  field  ,  i.e.,  ^  • 

The  rigid  rough  boundary  condition  was  perturbed  and 
statistically  averaged.  The  result  was  then  subtracted  from  the 
rigid  boundary  condition.  The  result  is  the  following. 


At  V*®  /  wMve 


(22) 


random  roughness 


random  slopes 


To  continue  the  stochastic  boundary  condition  into  space  from  the 
boundary.  Green's  formula  in  the  form  of  the  following  equation 
was  derived. 
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(23) 


A 


7«  -is*  ] 


Kuo"^  r 


and  the 

'in 


epresented 

reflected 


the  same  field  as 
field  ,  i  .  e .  , 


the  sum  of 


incident 


field 


♦(*=  (25) 

<t  =  JJexpC.(-4S-tit»3>J- 

•  [ 


Where 

*4,  *  ^  ^  i  -  ^  CXj*^y  . 


To  the  first  order  in  boundary  characteristics,  ^  will  be 
shown  later  to  be  the  sum  of  the  incident  field  potential,  , 

and  the  specularly  reflected  field  potential.  Therefore,  f  is 
the  total  xjotential  field  above  a  flat  boundary  located  at  the 
mean  rough  boundary.  This  identification  is  only  true  to  the 
first  order.  To  the  second  order,  the  specularly  reflected 
component  contains  roughness  def^ondent  terms  which  estimate  the 
energy  loss  from  specular  direction  to  all  scattering  directions. 

The  stochastic  field  potential,  ,  will  be  identified  as 
that  of  randomly  scattered  field  generated  by  the  boundary 
roughness  and  needs  statistical  analysis  for  its  physical 
implications. 
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Tho  form  of  ^  is  the  expanded  form  of  equation  (21)  in 
terms  of  roughness  height,  h,  and  slope,  s.  This 
particular  way  of  expansion  is  required  by  the  rigid 
boundary  condition.  Substituting  ^r  rough 

boundc'iry  condition,  followed  by  perturbation  in  term  of  the 
roughness  characteristics,  the  following  results  were  obtained; 


*  S  < A- ^  ^ X (27) 

K<«A„(*,>)=  jy  e*pf-^C4--S,V2'l' 

X' 

5  4  flo, « , JJ  ex  p  t- .•  i , )  •  5 ' 1  • 

(29) 


Therefore 


is  given  by 


=  exp£i  (■*,  *  +  *r3-)3 
^  J/exp[.-  tt-x + JJxxp  1- 
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(30) 


to  the  first  order  in  h  and  s.  The  previous  results  were 
modified  to  reflect  the  changes  in  the  incident  field  geometry, 
time  dependence,  and  the  change  in  the  boundary  conditin  to  rigid 
boundary . 

Accordingly,  from  equations  (24),  (25),  and  (30)  the 

following  expressions  can  be  obtained. 


f  s  ex)»£<«,  »-ilr|)J+  «pCi(-i,-a+Tfcr  J-)3 


(31) 


4-4 


=  JJexp[i  + 


(32) 
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Substitution  of  equation  (31)  into  (22)  results  in  the  expression 


2±! 


a  [-iVA-s 


By  utilizing  equation  (13),  equation  (23)  can  be  recast 


form: 


f  J|  ej<p  C'Ci  ’S  • 


(33) 


in  the 


(34) 
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Substitution  of  relationship  (33)  in  the  above  equation  will 
yield  the  same  result  as  in  equation  (32). 

This,  proves  the  equivalence  of  two  perturbation  methods  to 
the  first  order  in  the  random  roughness.  However,  the  total 
equivalence  requires  the  analysis  to  at  least  second  order  in  the 
statistical  properties,  e.g.,  random  roughness  and  random  slopes. 

4.  Conclusions 


Although  two  perturbation  scattering  theories  originally 

1  7 

formulated  by  Marsh  and  Bass  appear  to  have  the  same 
mathematical  approximation,  there  is  a  basic  difference. 

In  the  Bass's  method,  the  potential  field  is  perturbed  at  the 
start  by  utilizing  power  or  Taylor  series.  The  approximate 
solution  is  then  substituted  into  the  approximate  boundary 
conditions  perturbed  with  respect  to  the  random  roughness.  The 
perturbed  field  at  the  boundary  that  satisfies  the  approximate 
boundary  conditions  is  continued  into  space  by  either  using 
Green's  formula  (or  Helmholtz  integral)  or  using  its 
Quasi-Fourier  equivalent  form  of  the  quasi-plane  wave  expansion. 

In  the  Marsh's  method,  the  exact  solution  in  the  form  of 
quasi-plane  wave  expansion  is  substituted  into  the  appropriate 
exact  boundary  conditions.  The  scattered  field  at  the  boundary 
is  then  solved  by  utilizing  perturbations  in  terms  of  the 
boundary  characteristics  now  apparent  from  the  boundary 
conditions. 

The  Marsh's  method  is  preferred  because  perturbations  are  not 
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required  at  the  beginning  when  the  required  types  of 
perturbations  are  not  apparent. 
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